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TECHNICAL NOTE 3447 


ANALYTIC DETERMINATION OF THE DKCHARGE COEFFICIENTS OF FLOW NOZZLES 

By Frederick S. Simmons 


SUMMARY 

Integration of the velocity profile at the throat of a flow nozzle 
yields the discharge coefficient as a function of the ratio of houndary- 
layer thickness to the nozzle diameter. This ratio is obtained as a so- 
lution of the approximate momentum equation for the boundary layer. The 
resulting expression for the discharge coefficient is then a function of 
the Reynolds number based on nozzle diameter of the geometry of the 
nozzle. Good agreement is shown between this expression and published 
experimental data on flow nozzles for Remolds numbers between 10^ and 
106. 


INTRODUCTION 

When an ideal or frictionless fliiid passes through a nozzle^ the 
flow rate is a function only of the pressure drop, fluid properties, and 
nozzle geometry. For the measurement of the flow rate of an actual fluid, 
this functional relation must be modified to include the effects of fric- 
tion. This is usually done by introducing a "discharge coefficient", 
which is defined as the ratio of the actual mass-flow rate to the ideal 
mass-flow rate. Evaluation of this coefficient has been an experimental 
problem, and the results have usually been presented in the form of an 
empirical curve showing the discharge coefficient as a function of the 
Reynolds number . 

In rounded-approach nozzles with discharge coefficients close to 
unity, the frictional effects are concentrated in the boundary layer. A 
method of obtaining an analytical relation among the discharge coefficient, 
Reynolds number, and "Uie nozzle geometry by utilization of elementary 
boundary-layer theory is presented herein. 


ANALYSIS 


The actual rate of mass flow of a fluid through a nozzle is 
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^ ( 1 ) 

The integration will he taken over the eross-sectional area at the plane 
through the point of downstream pressure measurement, this plane herein- 
after being termed the "nozzle exit" . All symbols are defined in appendix 
A} the notation Is Illustrated In figure 1. 

The ideal mass-flov rate that would exist if the flow were completely 
Isentropic is 


% ” P%n^ ^2) 

The discharge coefficient for the case where the density is constant 
through the cross section is defined by 


C 




(3) 


In order to evaluate this integral, the velocity profile 

known over the cross-sectional area. If it is assumed that the boundary- 
layer thickness is small compared with the radius of the nozzle, then the 
exact shape of the velocity profile in the boundary layer need not be 
known. For integration of equation (3) , any of several functions could 
be used as approximations. The Blaslus solution (ref. 2) is the best 
known of these. However, for mathematical convenience in appendix B, 
which deals with the flow of gases with heat transfer, the function 


^ * tanh (ay) 

was selected and will therefore be used throughout this analysis. It 
will be shown later to be a sufficiently close approximation to the ve- 
locity profile of a laminar boundary layer for purposes of this 
investigation. 

A boundary- layer thickness B is defined by 

y at B for u/uq = 1 - e where 8 <<1 


so that 


^ M tanh ^b^^ where b » tanh“^(l-6) 


(4) 
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Upon substitution of equation (4), written for the nozzle exit, 
equation (3) may be integrated; however, the same result is more easily 
obtained by introduction of the displacement boundary- layer thiclmess 5*, 
which is defined by the relation 



(5) 


Since it is assumed that 8* << r, the boundary- layer build-up in a 
nozzle may be treated as a two-dimensional flow. The boundary layer 
produces a blocking effect equivalent to the reduction in area caused by 
the displacement thickness 8*, so that the effective exit area of a noz- 
zle with radius r is 


Ae,n “ 


m 3tr‘ 





( 6 ) 


The discharge coefficient is therefore given by the relation 


G * 1 



( 7 ) 


Substitution of equation (4) in equation (S) yields upon integration 


a I In 2 


( 8 ) 


Therefore, 


C 


2 In 2 ®n 
b r 


(9) 


The ratio 6u/r, of course, depends on the flow conditions and nozzle 
geometry. Consider now von Kanum’s momentum equation for the boundary 
layer for incompressible flow (ref. l) ; 



du^ 


u / u \ , 



dy a T 


( 10 ) 


where Uq is the velocity just outside the boundary layer, 
stress at the wall t is given by 


The shearing 
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Upon substituting equation (4) and performing the Integrations, 
equation (lO) becomes 

^ + 1 - iQ 2)^ “o ^ “ p(l - In 2) 


( 11 ) 


Equation (ll) may be Integrated after multiplying through by 2ug 

, „ In 2 

where n =s 2 ( 2 + ^ 


n-2 


6 


In 2 
2 


, with the result 
2 


2y^ 


P(1 - Ui 2) 




dx 


(12) 


4 


The solution of equation (12) for the case where Uq ** constant 
will be required later and is 



I 2 




'1 - In 2 

(13) 


i 




Expressed in terms of the displacement thlctoaess by substitution of 
equation (8), this solution is 

8f; _ 1.77 

The Blaslus solution for this case (ref. 2) gives 


^ = 1-75 
* 

The close agjreement Justifies the use of equation (4) to approximate the 
velocity profile. The two profiles are shown with another frequently used 
approximation in a dimensionless plot in figure 2. 

For the general solution of equation (l2) the function Uq(x) is 
required. Consider the convergent section of the nozzle shown in figure 
1. An approximation may be obtained for Uq(x) for this section by 
assuming the flow to be one-dimensional and incompressible. For 

0 < ^< |- and I’ « D, by continuity 


r 



NACA TN 3447 


5 


- 



- 2 sin 



(14) 


If this expression is used, equation (12) may he numerically inte- 
grated between the limits x/l ’ =0 and x/l ' = jt/2 giving 


I' 




(15) 


If the boundary layer is assumed to start at x/Z’ ■ jt/4 rather 
than at x/Z' »* 0 and equation (12) is integrated between the limits 
x/Z ' = jf/4 and x/Z’ = «/2, the value of the constant in equation (15) is 
not appreciably changed; this indicates that the flow very near the en- 
trance has little influence on the boundary- layer growth and thereby tends 
to justify the use of equation (14) . Moreover, if the linear relation 
Uq » Uq ^(2x/rtZ*) is arbitrarily assumed and substituted in equation (12), 

the value of the constant in equation (15) turns out to be 0.43, further 
indicating that accurate knowledge of the function Uq(x) is not important. 


The ratio &n/r for the entire nozzle is obtained in the following 
manner: The free-stream velocity Uq is assumed to vary according to 

equation (14) in the convergent section and is assumed to be constant and 
equal to Uq in the straight section. The equivalent straight section 

of length x', which produces the same boundary- layer thickness 6 as 
that produced by the convergent section, may be obtained by con^iarlng 
equation (15) with equation (13) written for x', with the result 

X' « 0.23 Z‘ 

Defining x = x' +2 and writing equation (13) for x yield 


5n 

Z' 


b A ~a/o-23 + -L 

A/l - In 2^ 2* 



or 
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Substitution of tMs relation in equation (3) gives for the discharge 
coefficient 


C * 





(17) 


DISCIBSION 

Discharge coefficients for the flow nozzles described in reference 3 
were calculated from equation (17). These values are shown, in figure 3 
compared with the e^erimentally determined values reported in reference 
3. In the range 10 < Rejj < 10°, the theoretical curve is well within 
the region of probable error of the experimental curve; this probable 
error in discharge coefficient ranges from +0.002 at Rej^ » 106 to +0.01 
at Rej) *10^. Theoretical curves obtained in reference 4 frcm a deter- 
mination of friction factors are also shown in figure 3. 

It is to be expected that equation (17) is valid for only a limited 
range of Reynolds numbers. For nozzles of the geometry considered in this 
analysis (Z/D or 2’/D, or both, of the order of unity) at low Reynolds 
numbers (Rejj < 10^), equation (lO) becomes less valid since it presupposes 
the boundary- layer thickness to be small compared with other lengths 
involved. At high Reynolds numbers (R^ > 10®) the boundary layer becomes 
turbulent, and equation (lO) again becomes less valid. 

Equation (lO), as written, is strictly valid only for incompressible 
fluids; however, the terms that would be added to equation (lo) to express 
the effects of compressibility are numericeLlly small compared with the 
terms already in that equation. Moreover, these terms would affect only 
the boundary- layer development in the convergent section of the nozzle. 
Hence, compressibility should have only a small effect on the discharge 
coefficient, and equation (17) may be considered valid for gases as well 
as liquids. 

Equation (l7) , of course, becomes inapplicable when extremes in noz- 
zle geometries are approached: When Z/d and I'/d become very large, 
the flow is fully developed; when Z/D and I’/D are very small, the 
flow approaches that throixgh a sharp-edged orifice. 

For maximum accuracy in a fluid-flow measurement with a flow nozzle, 
a calibration of course would be required; however, equation (l7) may be 
of value in extending a calibration curve that contains experimental dis- 
charge coefficients through only a small range of Reynolds numbers. 
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The analysis suggests an approach to the problem of determining 
changes in the discharge coefficients of flov nozzles for use in measuring 
the flow of gases when there is heat transfer between the gases nnA the 
nozzle surfaces. Thus, equation (3) could be written for the case where 
the density is variable through the cross section, by assuming a temper- 
ature profile as well as a velocity profile. A onide assumption con- 
cerning the effect of heat transfer on the boundary- layer growth would 
appear to be adequate to show the order of magnitude of -tlie changes to be 
expected in the discharge coefficient under these conditions. The anal- 
ysis for this case is presented in appendix Bj however, the validity of 
the result has not been established by experiment. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, February 14, 1955 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report: 
cross-sectional, area of nozzle 

tanh~^(l-e) b 
6 " 6 

t anh ~3-( l-e) 

discharge coefficient 

diameter 

length of straight section of nozzle 
axial length of convergent section of nozzle 

static pressure 
gas constant 
Reynolds number 
radius 

static teiT5)erature 
velocity 
mass-flow rate 

distance along nozzle surface 

equivalent straight length for convergent section of nozzle 
x« + 2 

distance perpendicular to surface of nozzle 


It 
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& Taoundary-layer thickness 
8* displacement houiidary-layer thickness 
6 a small fraction 

p viscosity 

p density 

T shearing stress at nozzle surface 

Subscripts: 
a actual 

c at end of convergent section of nozzle 
e effective 

h with heat transfer 
i ideal 

n nozzle exit 

o outside of houndary layer 

w nozzle vail 
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APPENDIX B 

DISCHARGE COEFFICIENTS FOR A GAS AT A TEMPEElATURE DIFFEEIEHT 
PROM THAT AT NOZZLE WALL 

The actual rate of mass flow of a fluid through a nozzle is 


We 


For a gas that obeys the law 


equation (Bl) may he written 


/pu dA 


pRt 


<B1) 




dA 


(B2) 


The ideal mass -flow rate for this case is 


- 1 ^ ^ 


o,n 

and the discharge coefficient is defined hy 

tr 


C«-f 

%n 


(B3) 


For the evaluation of this integral both a velocity and a temperature 
profile will be assumed, with the forms 


» tanh ( ay) 


(B4) 


where y « 5 for u/u^ ■ 1 - s and 

t - t^ 


■^o ■ "^w 


tanh (ay) 


(B5) 


t - t. 


where y »« 0 for t « t^ and y » B for 


1 - s. 


r~. 
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The assuiiiption of the same boundary- layer thickness 8 and the same 
shape for both the velocity and temperature profiles is essentially the 
Reynolds analogy. In each case the constant a, obtained by insertion of 
the boundary conditions^ has the value 


tanh"^ (l-e) _ b 
5 S’ 


Upon substitution of eq_uations (B4) and (B5) y eq_uation (B3) may be 
Integrated^ again^ hovever, the same result is more easily obtained by 
en^Jloyment of a displacement boundary- thickness. For this case 
the displacement thickness 8* vill be defined by 



(B6) 


Upon substitution of eq^uations (B4) and (B5) and integration^ 

to 

In 2 - In :r2 
^*8 ^ 

®h“ b 1: — 


2 - 




By analogy to the case of no heat transfer. 


C 



(B7) 


Therefore, 


C 


^ ^ lD.(io,n/'^) 

2 In 2 °n 1 TrTs 

-b— - 


(B8) 


If the assumption is made that the boundary-layer growth is not 
appreciably altered by the heat transfer and if the value of 8jj/r for 
the case of no heat transfer (eq.. (16)) is inserted in equation (B8), 
the discharge coefficient is given by 


4 In 2 




In 2 




€ 


1 II 
4 D 



C =* 1 - 


(B9) 
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where . 2 - 

The following table shows the variation of f(tQ^^/t^) with tp^^yt^: 


< 

o 

-p 


0.1 

4.92 

.25 

1.72 

.5 

1.33 

1 

1.00 

2 

.72 

4 

.50 

10 

.29 
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Figure 2. - Conparlson of various velocity profiles. 
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